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nonuniform  fields.  The  theory  is  based  on  the  use  of  similarity  and  related  arguments  to  bypass 
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IONIZING  POTENTIAL  WAVES  IN 
PREIONIZED  GASES 


1.  INTRODUCTION 

As  is  well  known,  the  application  of  a  high  potential  across  a  pair  of  electrodes  embedded  in  a 
weakly  ionized  gas  causes  waves  of  ionization  to  traverse  the  gas.  Ionization  occurs  primarily  at  the 
advancing  wave  head  where  a  high  electric  field  heats  existing  electrons  to  temperatures  sufficient  for 
electron-impact  ionization.  Impact  ionization  rapidly  raises  the  electron  density,  which  is  termed  an 
electron  avalanche,  until  the  resulting  wave  conductivity  becomes  so  large  that  space  charge  nullifies 
the  local  electric  field.  The  wave  thus  converts  field  energy  immediately  ahead  of  the  wave  into  gas 
ionization  and  excitation  energy  within  the  wave  body. 

The  advance  of  a  wave  is  governed  by  the  interplay  between  the  local  electric  field  £  and  the  wave 
conductivity  <r.  £  determines  the  growth  of  a-  due  to  electron-impact  ionization  while  <r  determines 
the  evolution  of  £  via  Maxwell’s  equations.  An  important  parameter  affecting  wave  propagation  is  the 
initial  conductivity,  a-,,  which  exists  immediately  ahead  of  the  wave.  A  number  of  mechanisms  have 
been  proposed  for  generating  this  conductivity  in  gases  that  are  initially  un-ionized.  Among  these 
mechanisms  are:  electron  drift,14  electron  pressure."  electron  runaway,8’ ?  gas  photoionization.8'13  and 
the  photoelectric  effect.14  The  proponents  of  these  different  mechanisms  all  claim  agreement  with 
experiment,  depending  upon  the  interpretation  of  the  experimental  results.15’ 16 

In  this  paper  we  circumvent  the  controversy  regarding  the  generation  of  a-,  by  considering  waves 
which  propagate  in  preionized  gases:  here  <r,  is  specified  as  an  initial  condition.  Such  an  analysis  applies 
to  all  but  the  first  ionizing  wave  propagating  through  a  gas.  The  extension  to  propagation  in  un-ionized 
gases  is  also  discussed. 

We  restrict  the  problem  to  nonuniform  fields  as  illustrated  in  Fig.  1.  At  time  zero  an  impulse  vol¬ 
tage,  «V  is  applied  to  a  sharply  contoured  electrode.  This  voltage  causes  an  ionizing  potential  wave 
(IPW)  to  emerge  and  propagate  toward  a  distant  ground.  The  fields  far  ahead  of  the  wave  are  negligi¬ 
bly  small.  Our  objective  is  to  determine  the  propagation  velocity,  K  of  the  IPW  as  a  function  of  <bn  and 
o-,. 


Previous  investigators  have  concentrated  on  determining  the  dependence  of  V  on  the  initial  elec¬ 
tron  density,  ntl,  which  is  directly  related  to  cr(.  Cravath  and  Loeb1  predicted  a  weak  logarithmic 
dependence  on  nei.  Schonland18  predicted  thet  V  scales  as  ne)/}.  Suzuki14  claimed,  with  some  experi¬ 
mental  justification,  that  V  scales  as  nj/2.  None  of  these  theories  indicates  the  dependence  of  V  on 
electrode  potential,  <£0-  Furthermore,  all  of  the  solutions  presented  for  V  depend  upon  one  or  more 
un-  known  variables.  These  variables,  which  include  the  wavefront  thickness  and  the  electron 
avalanche  rate,  may  themselves  depend  upon  nei  [see  also  the  criticisms  of  Winn20  and  Uman2ll. 

The  present  paper  attempts  to  remedy  these  shortcomings.  Two  objectives  have  been  imposed  on 
this  attempt.  The  first  is  to  obtain  solutions  which  depend  only  on  externally  specified  variables  such  as 
<&o  and  <rr  Explicit  dependence  on  the  'lectron  avalanche  rate,  for  example,  is  considered  unacceptable 
because  this  parameter  can  vary  by  orders  of  magnitude  and  a  unique  value  cannot  readily  be  chosen. 
The  second  objective  is  to  emphasize  the  physical  aspects  of  the  problem,  particularly  with  regard  to 
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obtaining  scaling  relationships.  As  a  result,  several  solutions  are  actually  presented.  The  first  solution 
is  essentially  a  dimensional  argument  which  minimizes  the  mathematics.  Subsequent  solutions  indicate 
limitations  and  modifications  which  arise  in  certain  parameter  regimes  and  applications. 


2.  A  SIMPLE  SCALING  MODEL 


Consider  the  basic  problem.  An  electtode  is  held  at  constant  voltage.  <j>{).  A  long  filamentary 
IPW  emanates  from  the  electrode  and  propagates  with  velocity  V  through  a  homogeneous,  weakly  ion¬ 
ized  gas  of  conductivity  cr,.  We  assume  that  f’is  determined  solely  by  and  <r,.  Dimensional  analysis 
then  demands  that 

P  -  So  <t>n  \ncrJE,  (1) 

where  £,(<£, >.  <r,)  is  a  characteristic  field  strength  and  go  is  a  numerical  constant. 

To  determine  £,(<t>o-  o-,),  we  note  that  two  principal  rates  govern  the  behavior  of  IPWs:  47r<r 
which  controls  the  evolution  of  field  £;  and  the  electron  avalanche  rate.  5.  which  controls  the  growth  of 
conductivity  <r.  S  is  principally  determined  by  £and,  to  a  lesser  extent,  by  the  degree  of  ionization  and 
excitation.  For  simplicity  £  is  assumed  to  be  a  known  function  of  £and  <r,  where  wave  conductivity  a 
is  a  measure  of  degree  of  ionization. 

We  can  thus  identify  two  characteristic  rates:  4jt<7,  and  £(£,.  cr,).  We  presume  that  these  two 
rates  are  related  and  set  them  equal  to  within  a  numerical  coefficient,  g f; 

£(£,.  cr,)  *  gi  4rrcr,.  (2) 

This  equation  determines  £,  for  a  given  initial  conductivity  a,.  Equation  (1)  then  determines  wave 
speed  V.  Estimates  of  £(),  ,?j  are  presented  later. 

To  understand  the  physical  origin  of  these  relationships,  imagine  doubling  all  physical  dimensions 
of  the  IPW  shown  in  Fig.  1  while  simultaneously  doubling  the  electrode  voltage  d>0;  i.e..  all  characteris¬ 
tic  lengths  I,  become  21,  while  <60  becomes  2<60.  This  doubling  redistributes  but  otherwise  does  not  alter 
the  local  variables  such  as  <r.  S,  and  £  (which  scales  as  (£<//,>•  Since  the  characteristic  rates  are  unal¬ 
tered  while  the  characteristic  lengths  double,  wave  speed  V  which  scales  as  SI,  also  doubles.  We  thus 
conclude  that  V  scales  linearly  with  <t>{)  and  that  £f  in  Eq.  (1)  is  a  function  only  of  cr,.  The  important 
point  is  that  in  the  absence  of  geometrical  constraints  on  length  scales  the  initial  conductivity  <r, 
determines  the  field  E,  while  the  applied  potential  <60  controls  the  range  over  which  E,  is  distributed. 

Equations  (1)  and  (2)  form  the  foundation  of  the  present  theory.  What  follows  is  an  examination 
of  more  subtle  issues.  These  include:  (i)  derive  a  solution  based  on  a  closed  set  of  wave  equations; 
(ii)  identify  and  incorporate  other  wave  parameters  of  interest;  and  (iii)  determine  when  the  model  fails 
or  needs  modification. 


3.  WAVE  EQUATIONS 


We  now  present  a  more  formal  derivation  of  Eqs.  (1)  and  (2)  and  assess  some  of  their  limita¬ 
tions.  We  express  Maxwell's  equations  by 


£-  -V  <6 


(3) 


and 


V  • 


+  4  irtrE 


0 


(4) 
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where  <t>  is  the  electrostatic  potential.  Ignoring  electron  diffusion  and  ignoring  the  massive  and  essen¬ 
tially  immobile  ions,  the  wave  conductivity  is  given  by 

<t  —  enefi  <  5 ) 

where  e  is  the  electron  charge,  ne  the  electron  density,  and  the  electron  mobility.  The  electron  den¬ 
sity  evolves  according  to  the  continuity  equation 

~  ne  +  V  ■  (ne  u)  »  S  ne  (6) 

where  u  is  the  electron  fluid  velocity  and  where  the  electron  avalanche  rate  S  equals  the  net  electron 
production  rate  due  to  collisions.  In  the  absence  of  electron  diffusion,  u  equals  the  electron  drift  velo¬ 
city, 

u  —  —fi  £.  (7) 


Equations  (3)-(7)  are  dosed  by  specifying  the  avalanche  rate  S  and  mobility  n  in  terms  of  the 
field  E.  the  conductivity  cr,  and  the  gas  density  .V.  On  the  time  scales  of  interest  we  assume  that  .V 
does  not  change  so  that  5  and  n  are  functions  only  of  £  and  cr  : 

S-S(£.o-)  (8) 

H  «  n  ( £.  cr ) .  <  9 ) 

For  a  very  weakly  ionized  wave  the  dependence  on  cr  would  also  vanish;  at  higher  degrees  of  ioniza¬ 
tion,  processes  such  as  electron-ion  recombination  and  super-elastic  collisions  with  excited  molecules 
affect  the  electron  dynamics,  and  thus  affect  the  electron  transport  coefficients. 

Equations  (3M9)  plus  the  appropriate  boundary  conditions  determine  the  evolution  of  ionizing 
potential  waves.  To  obtain  stationary  solutions  in  a  wave  frame  traveling  with  constant  velocity  V  with 
respect  to  the  laboratory  (gas)  frame,  we  use  the  traveling-wave  condition, 

| - V-V  (10) 

at  ~  ~ 

to  rewrite  Eqs.  (4)  and  (5)-(7)  as 

V  •  ([_K-  Vl  £  -  4ir<r  £)  -  0  (11) 

and 

V  •  ([_F+  /u£]  crlix)  ”  —  Srr/ju.  (12) 


To  obtain  scaling  relationships,  we  look  for  self-similar  solutions  of  the  form 

q(x,  r)  *•  (13) 

where  q  is  any  variable,  q,  is  a  characteristic  value,  and 

£*»<£—  Vt)i !,  (14) 

with  I,  being  a  characteristic  length.  The  functions  /?(f)  are  interrelated  through  the  atomic-physics 
relationships  (8)  and  (9)  and  through  a  dependence  on" cr,  which  characterizes  the  initial  state  of  the 
gas.  The  functions  do  not  depend  on  electrode  potential  <£0.  Note  that  Eqs.  (13)  and  (14)  explicitly 
incorporate  the  traveling-wave  condition  (10). 


Using  the  transformation 


4 


we  reduce  Eqs.  (3),  (11),  and  (12)  to 


(16) 


E,  -  C<V/„ 


b 

V 

a 

d/a 

,  d/4 

1 ,  4jtct, 

‘ 

& 

K  if 

a 

F+m£ 

fa 

fa  fs 

if 

I'S, 

4 

4 

where  C.  is  a  dimensionless  constant  and  where  the  characteristic  avalanche  rate  is  taken  to  be 

S,  =  S(£„  a-,). 


(17) 


(18) 


(19) 


Equations  ( 16)-(  19)  preserve  self-similarity  provided  the  terms  in  square  brackets  are  constant. 
This  is  valid  provided  Eqs.  (1)  and  (2)  are  satisfied,  and  provided  the  electron  drift  velocity.  —  fiE,  is 
negligibly  small.  Since  the  electrostatic  Eq.  (3)  is  justified  only  when  the  wave  speed  is  small  compared 
with  the  speed  of  light,  c.  we  expect  Eqs.  (1)  and  (2)  to  be  valid  only  for 

\(iE, I  <  V  <  c.  (20) 

The  self-similar  analysis  and  Eqs.  (1)  and  (2)  become  invalid  at  low  wave  speeds.  This  failure 
arises  because  the  electrons,  which  are  the  dynamic  particles  controlling  wave  conductivity,  are  not  sta¬ 
tionary  in  the  laboratory  frame  but  are  drifting  with  velocity  ~/iE.  The  impact  of  this_drift  can  be 
qualitatively  understood  by  ignoring  radial  drift  and  assuming  constant  axial  drift,  -nE,.  Equations 
(17M19)  then  require  that 

5(£,,<r,)  —  ,C|4ircri(l  4-^i£,/F).  (21a) 

Or  equivalently, 

g,  E,  _ 

S(£,,<t,) - —  /x£,  —  g |  4t rcr,,  (21b) 

So  ®o 

replacing  Eq.  (2).  The  latter  equation  demonstrates  that  £,  is  in  general  determined  by  (60  as  well  as  by 
(T,.  Furthermore,  positive  waves  for  which  <60  >  0  require  a  larger  |£,l  than  negative  waves  for  which 
<bn  <  0.  According  to  Eq.  (1),  positive  waves  should  therefore  propagate  slower.  Indeed.  Eq.  (21a) 
rigorously  restricts  the  negative-wave  speed  to  V  >  |/*£,|,  which  reflects  the  fact  that  a  wave  cannot 
propagate  more  slowly  than  its  constituent  particles.  A  similar  restriction  applies  to  positive  waves  only 
if  S(£,,<r,)  increases  no  faster  than  linearly  with  £,. 

This  section  has  provided  a  formal  basis  for  the  IPW  velocity  solution,  V(<t>0,  <r,).  presented  ear¬ 
lier.  The  range  of  applicability  has  also  been  discussed.  The  next  section  takes  yet  a  different 
approach,  discusses  additional  limitations,  and  provides  estimates  for  the  coefficients  g0,  gt.  and  C.. 

4.  TRANSMISSION-LINE  MODEL 

To  simplify  the  problem  further,  we  use  a  transmission-line  analysis  to  approximate  Maxwell’s 
equations.  The  electric  field  along  the  axis  of  propagation  is  then  given  by 

£  -  -  -  4-  (L-7)  <22) 

b:  bt  ' 

where  the  last  term  allows  for  inductive  corrections.  Ohm's  law  relates  the  current  I  to  the  wave  con¬ 
ductance  1  via 


/-!£  (23) 

where  I  equals  a  radial  integral  of  wave  conductivity,  cr.  Equation  (12)  yields 

-■|-  [(»' +M£>  X//*l- SX/M  (24) 

o: 

where  £  is  a  radial  average  of  5. 

Charge  conservation  produces 

(CM  +  |^  -  0  (25) 

3/  ■  d: 

where  C.<t>  is  the  space  charge  distributed  per  unit  length  along  the  wave.  The  distributed  capacitance 
C.  and  inductance  L.  satisfy 

(L.  C.rl/:  <  c.  (26) 


We  assume  for  convenience  that  (£.  C.)  is  constant  and  use  the  traveling-wave  condition  (10)  to 
obtain 


£--  32.  (i  _  F-Z..C) 
3r 

/  -  vc.  6  -  I  £ 


(27) 

(28) 


and 


f:  In  (M£) 


CAV  +ixE) 


f  +  m  Ed,  ,,  S 


(29) 


1(1  -  K-L.C) 

Equations  (26>-(29)  demonstrate  that  inductive  effects  serve  principally  to  limit  the  wave  speed  to 

V  <  (£.C.)-I/*  <  c.  (30) 

Henceforth  we  generally  ignore  inductive  effects  except  to  impose  constraint  (30). 


Let  us  use  Eqs.  (24)  and  (29)  to  describe  the  qualitative  structure  of  IPWs.  Well  within  the  wave 
the  conductance  X  approaches  its  maximum  value.  X;  .  while  the  field  £  approaches  a  low  value.  £,.  £, 
is  a  maintenance  field  defined  by  the  condition 

5<£r,  <rf)  =  0  (31) 

where  a- ,  is  the  final  conductivity  attained  in  the  wave  body.  If  <t;  were  zero.  £,  would  equal  the 
breakdown  field,  Eb:  for  large  cr,,  E,  can  be  well  below  Eb  due  to  the  ionization  and  electron  heating 
provided  by  the  excited  molecules. 


Approaching  the  head  of  the  wave,  £  and  S  begin  to  rise  while  X  gradually  falls.  As  S  becomes 
large,  the  fall  of  X  accelerates.  As  a  result,  £  and  5  attain  peak  values.  £,  and  S,,  respectively,  at  the 
wave  tip.  Beyond  the  tip  £  and  S  approach  zero.  In  this  region  the  distributed  space  charge,  C-d>.  and 
capacitance.  C,  also  approach  zero. 


The  transmission-line  equations  are  incomplete  since  the  wave  radius.  r„.,  and  the  axial  depen¬ 
dence  of  C.  are  unknown.  We  thus  supplement  these  equations  by  modeling  the  IPW  as  a  long, 
space-charged  filament.  Electrostatics  then  dictates  that  axial  variations  in  the  wave  head  are  character¬ 
ized  by  a  common  length  scale, 

/,  -  C.  6./E,  (32) 


6 


where  C.  and  0,  characterize  the  capacitance  and  potential,  respectively,  of  the  wave  head.  Typically. 
C.  =  0.2  in  Gaussian  units.  Note  that  a  common  length  scale  can  be  imposed  only  if  condition  (20)  is 
satisfied.  See  Appendix  A  for  details. 


We  express  the  tip  conductance  by 


I,  -  g‘  <r,irr; 

(33) 

and  set  the  wave  radius  to 

(34) 

The  strong  dependence  of  the  avalanche  rate  S  and  weak  dependence  of  the  mobility  fj. 
gests  that  g'  =  1  and  g''  <  1. 

on  field  E  sug- 

Inserting  these  assumptions  into  Eqs.  (24)  and  (29)  and  imposing  constraint  (20) 
lined  in  Appendix  A, 

yieids.  as  out- 

f  -  0  f  -)7T<T  ,/  E, 

(35) 

and 

S(E,.  <r,)  =  S,  =*  gj  417(7, 

(36) 

where 

go  “  g'(g"C.)2/4C:r 

(37) 

and 

.?!  “  SiJC.. 

(38) 

Here  the  capacitance  at  the  wave  tip  satisfies  C-.  «  C.. 


Equation  (35)  is  equivalent  to  Eq.  (U  provided  the  tip  potential  0,  equals  the  electrode  potential 
0O-  This  equivalence  is  justified,  to  order  C.  «  1.  provided  the  resistive  voltage  drop  within  the  wave 
body  is  small.  This  proviso  is  also  required  to  justify  the  assumption  of  constant  velocity  V.  \  generally 
sufficient  condition  is  that  the  IPW  length  satisfies 

I  «  /maK  -  0o/£>  (39) 

where  /ma,  is  the  nominal  maximum  propagation  distance.  We  point  out  that  IPWs  can  actually  pro¬ 
pagate,  with  diminishing  velocity,  for  distances  further  than  /ma,  due  to  the  space  charge  stored  in  the 
wave  body.  This  charge  and  the  capacitive  structure  of  IPWs  permits  propagation  for  several  "RC" 
times,  suggesting  that  distances  up  to  3/mas  may  be  possible.  This  feature  additionally  explains  the  abil¬ 
ity  of  IPWs  to  continue  propagating  even  if  the  electrode  voltage  is  suddenly  terminated." 

The  transmission-line  analysis  and  electrostatic  model  have  reproduced  the  original  solutions 
given  by  Eqs.  (1)  and  (2).  We  now  use  these  models  to  estimate  the  coefficients  g0  and 

Equation  (35)  has  a  simple  physical  basis.  An  IPW  propagates  by  transforming  conduction 
current  in  the  wave  body  into  displacement  current  at  the  wave  tip.  This  transformation  produces  high 
space-charge  fields  at  the  tip  which  generate  high  conductivity  in  the  body.  Conduction  currents  are 
negligible  at  the  tip  only  if  the  electric  field  rises  to  its  peak  value  E,  in  a  short  distance 

I,  «  K/4rcr,.  (40) 

Equation  (32)  converts  this  condition  into 

V'»  C.<t>,  4w(r,/f,.  (41) 
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The  strength  of  this  inequality  is  limited  by  the  fact  that  conduction  currents  dominate  inside  the 
wave.  Hence.  Eqs.  (35),  OS),  and  (41)  suggest  that 

j?o  —  1 0  C.  =  2  (42) 

and 

.if,  =  10.  (43) 

To  illustrate  these  findings,  consider  very  weakly  preionized  air  of  density  .V.  The  reduced 
avalanche  rate.  S/N.  is  then  well  known23  ’-  as  a  function  of  reduced  field.  E/N,  as  shown  in  Fig.  2.  A 
corresponding  plot  of  V/<t>n  as  a  function  of  a, /.Vis  shown  in  Fig.  3.  The  latter  plot  is  based  on  Eqs. 
(1).  (2).  (42),  and  (43).  Also  plotted  in  Fig.  3  is  the  electron  drift  speed.  l/a£,i,  which  is  presumed 
small  compared  with  wave  speed  )'. 

Several  features  of  these  plots  are  noteworthy.  First,  as  <r,/.V  becomes  small  E,  approaches  Eh, 
where  the  breakdown  strength  EJ  S  =  100  Td  in  air.  [A  Townsend  <Td)  equals  10_I  V’olts-cm2  — 
0.33  Volts/cm-torr.]  As  a  result  K/«40  scales  nearly  linearly  with  crJS.  This  dependence  lessens  at 
high  <7,/.V.  At  <tJS  >  10-l°  cm3/sec.  F'/c6o  is  essentially  constant.  Here  S  approaches  which  is 
the  maximum  possible  ionization  rate  in  the  gas.  In  very  weakly  preionized  air.  Smux/.V  =  5  x  10“* 
cmVsec. 

The  restriction  on  avalanche  rate  S  indicates  that  relationship  (2)  is  valid  only  for 

O’,  <  O’ max  =  Smix/4irg,.  l44) 

Above  this  limit  conduction  currents  ahead  of  the  wave  redistribute  the  wave  space  charge  and  associ¬ 
ated  fields  before  the  conductivity  can  avalanche.  We  thus  conclude  that  unbounded  traveling-wave 
solutions  cease  to  exist  at  high  a  .  This  failure  is  closely  related  to  the  energy  arguments  of  Appendix 
B  which  show  that 

E/<bi)  <5x10'  cm/Volt-sec. 


A  different  failure  arises  at  low  o-,  where  E,  approaches  Eh.  Here  electron  attachment  generates 
large  numbers  of  negative  ions  and  ion  currents  which  can  short-circuit  the  electron  avalanche  process. 
We  express  the  avalanche  rate  as 

5  =  (a  —  tj)  fiE  (45) 

where  a^E  is  the  electron-impact  ionization  rate  and  t )fxE  is  the  net  electron-attachment  rate:  a  is 
termed  the  first  Townsend  coefficient  and  17  the  attachment  coefficient.  Ion  currents  can  be  ignored 
provided  the  production  of  negative  ions  due  to  electron  attachment  does  not  exceed  the  production  of 
electrons.  We  thus  require  that  the  peak  avalanche  rate  satisfies 

S.  >  r)(Fi,)fiEh  (46) 

and  hence  that 

o’,  >  o’ mm  =  V(Eb)(*Eb/4ngl  (47) 

where  the  breakdown  field  Eb  is  defined  by 

a(Eb)  -  (48) 

We  identify  <rmin  as  the  minimum  initial  (electron)  conductivity  required  to  augment  1PW  propagation. 
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)  N/S 


E/N  (Td) 


Fig.  2  —  Electron  impact  ionization  rate.  S.  as  a  function  of  field  strength.  E.  in  weakly  ionized  air  of  density  of  V 
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In  cold  weakly  preionized  air,  t}(£4)/.V=  1.5  x  10  19  cm2  and^£6  =  1.5  x  107  cm/sec.  Assum¬ 
ing  g  —  10.  conditions  (44)  and  (47)  reduce  to 

10-u  <  —  <  5  x  10~10  cmVsec.  (49) 

/V 

Because  u ,  scales  as  the  degree  of  preionization.  nJS\  the  lower  limit  is  generally  more  important  at 
high  .V  while  the  upper  limit  is  more  important  at  low  S. 

5.  RADIALLY  BOUNDED  IPWs 


In  the  preceding  theory  the  wave  radius  rw  was  determined  by  6o  and  cr,.  We  now  consider  the 
modifications  which  arise  if  r„.  is  externally  constrained,  either  by  the  presence  of  walls  or  by  the  initial 
profile  of  cr,.  We  term  these  waves  radially  bounded. 


Radial  constraints  have  two  main  effects.  The  first  is  that  the  internal  wave  structure  is  altered  so 
that  a  single  axial  length  no  longer  characterizes  I,  £.  and  C..  In  particular,  the  distributed  capacitance 
C.  (rather  than  potential  4>)  varies  slowly  throughout  the  wave  head  and  may  be  treated  as  constant. 
Setting  C  —  C.  and  evaluating  Eqs.  (29)  and  (33)  at  the  peak  field  £,  yields 

V  -  g:  (47TO-,  S,)l/:  rh  (50) 

where  rh  is  the  bounded-wave  radius  and 


Si 


g  l  —  V2  L.  C. 
4  C  1  +HEJV 


A  nearly  identical  equation  was  obtained  by  Suzuki.19 


1. 


(51) 


The  second  effect  is  that  £,  is  determined  not  by  the  initial  conductivity  but  by 

£,  -  g"C.<t>o/rb,  (52) 

provided  the  resistive  voltage  drop  in  the  wave  body  is  negligible.  This  equation  is  a  direct  extension  of 
Eqs.  (32)  and  (34)  and  closes  the  velocity  solution  of  Eq.  (50)  by  defining  the  peak  avalanche  rate. 

S,  =  S(£„  a-,).  (53) 


The  transition  from  the  unbounded  wave  solutions,  given  by  Eqs.  (1)  and  (2).  to  the  bounded 
solutions,  given  by  Eqs.  (50M53).  occurs  when  the  bounded  wave  radius  rb  becomes  less  than  the 
unbounded  radius.  This  transition  typically  occurs  in  air  when 

<60 /r„N  >  10"13  Volt-cnr.  (54) 

The  character  of  radially  bounded  IPWs  differs  from  that  of  unbounded  IPWs.  Bounded  IPWs  are 
predicted  to  behave  as  follows:  wave  speed  V  scales  nearly  linearly  with  cr,l/2:  V  exhibits  a  monotonic 
but  generally  nonlinear  dependence  on  <60;  v  exhibits  a  nonmonotonic  dependence  on  bounded-wave 
radius  rb  and  on  gas  density  N.  Moreover,  constraints  (44)  and  (47)  on  r,  are  no  longer  valid  since 
Eq.  (2)  does  not  apply.  In  particular,  bounded-wave  solutions  exist  even  at  high  <r,  »  tZma,-  The  rea¬ 
son  is  that  bounded  waves  are  constrained  to  be  filamentary  which  is  the  characteristic  required  for 
space-charge  field  enhancement;  i.e.,  the  absence  of  conductivity  outside  rb  ensures  that  space-charge 
fields  and  subsequent  avalanching  will  develop  inside  rb. 
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6.  EXTENSION  TO  <r,  —  0 

For  many  problems  of  interest  the  initial  conductivity  ahead  of  the  wave  is  well  below  the  value. 
ermm,  required  to  influence  IPW  propagation.  Waves  can  still  propagate,  however,  provided  the  wave 
itself  can  generate  the  necessary  seed  ionization  at  the  wave  tip. 

The  only  mechanism  explicitly  incorporated  in  the  theory  for  generating  this  seed  ionization  is 
electron  drift.  Here  the  directed  motion  of  the  electrons  due  to  the  local  electric  field  advances  the 
wave.  Equation  (29)  indicates  that  for  vanishing  tip  conductance,  wave  speed  Fapproaches 

F  =  -  ix  E,-  (55) 

This  is  positive  definite  only  for  negative  waves  where  V  •  E  <  0.  Electron  drift  can  thus  account  for 
slow  negative  waves  but  cannot  account  for  fast  (F  »  108  cm/sec)  negative  waves  or  for  positive 
waves. 

The  theory  can  easily  accommodate  other  seed  ionization  mechanisms  provided  the  tip  conduc¬ 
tivity  or  conductance  can  be  estimated.  The  only  formal  constraint  is  that  condition  (47)  be  satisfied. 
A  less  formal  constraint  is  that  the  seed  mechanism  account  for  positive  waves  as  well  as  negative 
waves  since  the  experimental  character  of  both  waves  is  so  similar.  This  last  constraint  excludes,  for 
example,  the  electron-runaway  mechanism.6 

Two  proposed  mechanisms  in  nonuniform  fields  are  gas  photoionization* 13  and  electron  pressure.' 
(The  photoelectric  effect14  is  presumed  irrelevant  because  the  ground  electrode  is  in  a  weak-field  region 
far  from  the  wave.)  For  positive  waves  the  electron-pressure  mechanism  requires  an  anomalous 
electron-heating  process  which  so  heats  electrons  in  the  positive  wavefront  that  the  electron  fluid  velo¬ 
city  reverses  direction  and  opposes  electron  drift.  Such  a  reversal  is  deemed  improbable  because  these 
electrons  would  then  rapidly  lose  energy  to  the  electric  field  rather  than  gain  it. 

A  major  criticism  of  the  photoionization  mechanism  is  that  the  relevant  atomic  physics  is  often 
unknown.14  26  Photoionization  nonetheless  remains  promising  because:  (i)  ample  experimental  evi¬ 
dence27  demonstrates  that  it  or  a  similar  process  occurs  in  electrical  discharges  and.  by  inference,  occurs 
in  IPWs;  and  (ii)  inclusion  of  the  experimentally  determined  photoionization  rate  in  the  electron  con¬ 
tinuity  equation  directly  leads  to  fast  positive  and  negative  waves.1213 

To  illustrate  the  role  of  seed  ionization  mechanisms  on  wave  propagation,  consider  gas  photoioni¬ 
zation.  Ionizing  photons  are  created  within  the  wave  with  an  efficiency  8  per  electron-impact  ionization 
event.  The  photon-production  process  determines  whether  8  depends  on  gas  density  .V.  field  parameter 
E/S.  or  degree  of  ionization.  Typically,  8  «  l.8  Such  a  low  efficiency  means  that  photoionization  is 
unimportant  in  the  wave  head  or  body. 

Photoionization  can  nonetheless  be  important  ahead  of  the  wave  provided  the  photon-absorption 
length,  Ka,  is  much  longer  than  the  characteristic  electron  transport  lengths.  Photoionization  then  gen¬ 
erates,  within  a  length  Ka  outside  the  wave  head,  an  initial  conductivity  of  <r,  =  8  <t ,  where  <t ,  is  the 
final  wave  conductivity.  This  seed  ionization  is  sufficient  to  sustain  wave  propagation  provided  the  gain 
in  conductivity,  crj/cr,,  due  to  electron-impact  ionization  in  the  wave  head  exceeds  S-1.  See  Ref.  28  for 
further  discussion. 

7.  ADDITIONAL  LIMITATIONS  AND  MODIFICATIONS 

Two  additional  parameters  which  can  affect  F  and  thereby  invalidate  or  modify  the  theory  are 
electrode  effects  and  a  finite  rise-time,  r,,  of  the  applied  voltage.  Only  if 

~r  <  <  UJ  F  (56) 
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can  the  voltage  rise-time  be  ignored.  If  this  condition  is  not  satisfied,  the  effective  electrode  voltage  is 
reduced  and  Eq.  ( l )  should  be  adjusted  accordingly.  Alternatively  if  <£0  rises  linearly  with  time,  wave 
speed  is  ultimately  given  from  Eqs.  (10)  and  (27)  by 


V 


E7 


~JT 


(I 


V'-L.C.) 


<57) 


where  the  asymptotic  field  strength  Ef  is  defined  by  Eq.  (31). 


Electrode  effects  assume  a  particular  importance  when  examining  the  differences  between  positive 
and  negative  waves.  Positive  waves  generate  an  electron  flux  that  flows  into  the  electrode,  whereas 
negative  waves  require  an  electron  flux  that  emanates  from  the  electrode.  Negative  waves  thus  require 
an  active  electron-emitting  surface  whereas  positive  waves  require  an  essentially  inactive,  electron¬ 
absorbing  surface.  Any  restriction  on  electron  current  flowing  in  or  out  of  the  electrode  results  in  resis¬ 
tive  sheaths,  which  reduce  the  voltage  appearing  at  the  wave  head  and  thereby  reduce  wave  speed  V. 
Electrode  sheath  restrictions  are  thus  usually  more  severe  for  negative  waves  than  for  positive  waves. 
Such  restrictions  and  differences  become  unimportant,  however,  when  the  applied  potential  is  large 
(typically.  <t>o  >  10  kV)  compared  with  the  sheath  potential. 


8.  DISCUSSION  AND  COMPARISON  WITH  EXPERIMENT 


Let  us  briefly  review  the  previous  theoretical  predictions.  We  have  identified  two  types  of  waves: 
radially  bounded  and  radially  unbounded.  Unbounded  means  that  the  wave  radius  is  smaller  than  any 
external  constraint.  The  velocity  of  unbounded  waves  is  predicted  to  scale  at  low  <rJN  as 

V  -  <60  <r  J -V  -  4>nnJ N2  (58) 

where  .V  is  the  gas  density  and  na  is  the  initial  electron  density  ahead  of  the  wave.  By  contrast,  the 
velocity  of  radially  bounded  waves  is  predicted  to  scale  as  (r}/:  and  to  depend  in  a  complicated  manner 
on  <a0.  V,  and  bounded-wave  radius  rb.  At  high  <f>0  and  low  ,V,  rb  this  dependence  simplifies  to 

V  -  (<r(.V)l/J  r„  -  ne\/2r„.  (59) 

Electrode  effects  and  other  processes  can  easily  modify  these  scalings. 


Corresponding  to  the  two  types  of  waves  are  two  classes  of  experiments:  JPWs  launched  in  low- 
pressure  glow  discharges  where  the  initial  conductivity  is  determined  by  the  glow-discharge  current:  and 
high-pressure  experiments  using  lasers  to  preionize  the  gas.  A  third  pertinent  category  is  the  long 
spark29  in  which  a  succession  of  IPWs  (streamers,  leaders,  and  return  strokes)  propagate  along  the  same 
path.  A  difficulty  with  the  latter  experiments  is  that  the  state  of  the  gas  between  successive  waves  is 
unknown,  and  hence  or,  is  indeterminate.  We  point  out  that  each  successive  wave  encounters  an  initial 
conductivity  much  higher  than  that  experienced  by  its  predecessor.  As  a  result  wave  speed  is  higher 
and  wave  radius  is  (slightly)  smaller.  Successive  IPWs  are  therefore  predicted  to  be  radially  unbounded 
provided  the  predecessor  is  (until  <t,  exceeds  crma,  as  defined  by  Eq.  (44)). 


We  focus  on  a  limited  subset  of  the  availa'  !e  experimental  data.  Consider  first  the  laser-guided 
discharges  of  Koopman  and  Saum30  and  of  Greig,  et  al.31  32  These  experiments  used  a  pulse  of  energy 
from  a  Nd:  glass  laser  or  a  C02  laser  to  preionize  meter-long  paths  in  atmospheric  air.  The  subsequent 
application  of  an  impulse  voltage  across  electrodes  at  each  end  induced  IPWs  to  propagate  within  the 
laser-formed  channel,  ultimately  leading  to  a  high-current  discharge.  The  initial  conductivity  was  con¬ 
trolled  by  varying  either  the  laser  energy  or  the  delay  time  between  firing  the  laser  and  applying  the 
impulse  voltage. 
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The  minimum  laser-generated  electron  conductivity.  o-mm,  required  to  guide  IPWs  was  consistent 
with  predictions  (47)  and  (49).  Analysis  of  the  experiments  indicated  that  thermal  detachment  of 
negative  ions  by  laser-heated  molecules,31-32  rather  than  photodetachment  by  IPW  photons,30  played  a 
key  role  in  determining  crmjn. 

The  experiments  of  Greig,  et  al.  attempted  to  determine  the  dependence  of  wave  speed  V  on 
initial  conductivity  <r,  and  electrode  voltage  <j>n-  Wave  speed  was  determined  by  measuring  the  arrival 
time  of  the  potential  wave  at  several  axial  positions  along  the  laser-designated  path  using  capacitive  vol¬ 
tage  probes.  Results  for  Fas  a  function  of  <60  are  shown  in  Fig.  4  for  an  initial  conductivity  of  cr,/ N  = 
!0-12  cm3/sec.  Not  only  is  the  plot  linear,  but  the  measured  value  of  F/<£0  =  5  x  !03  cm/Volt-sec 
agrees  to  within  experimental  uncertainty  (a  factor  of  2  or  more  for  a-JM)  with  that  predicted  by  Fig.  3. 
Such  agreement  is  remarkable  considering  the  order-of-magnitude  estimates  used  for  coefficients  and 
j?,.  The  dependence  of  V on  <r,  was  less  clear  due  to  a  large  scatter  in  the  conductivity  measurements. 
Fwas  shown,  however,  to  exhibit  a  strong  monotonic  dependence  on  cr,. 

Several  other  features  of  Greig’s  experiments  are  noteworthy.  First,  optical  diagnostics  indicated 
that  the  IPWs  were  radially  unbounded  since  the  wave  radius  (/•„.  <0.1  cm)  was  smaller  than  the 
radius  (0.2  to  2  cm)  of  the  laser-formed  channels.  Second,  the  guidance  afforded  by  the  laser  preioni¬ 
zation  was  strong;  IPWs  could  be  guided,  for  example,  even  along  laser  paths  that  were  perpendicular 
to  the  static  electric  field  lines.  This  plus  the  large  increase  in  guided  wave  speed  over  unguided  wave 
speed  again  supports  the  prediction  that  cr,  strongly  influences  wave  propagation.  Third,  wave  speed 
and  wave  potential  decreased  with  increasing  propagation  length.  I.  And  fourth,  the  maximum  propaga¬ 
tion  length  for  a  given  electrode  potential  </>n  corresponded  to  an  average  electric  field  l)  of 

£min  =  1.5  kV/cm.  (60) 

This  field  is  roughly  one  third  of  the  maintenance  field.  £,.  predicted  by  numerical  simulations  based 
on  the  detailed  air  chemistry  code.  CHMAIR.33  These  simulations  assumed  that  laser  preheating  raised 
the  channel  temperature  from  300°K  to  800°K  and  lowered  the  gas  density  accordingly.32  The  many 
approximations  employed  in  the  simulations  preclude  an  accurate  determination  of  £,.  The  results 
nonetheless  support  the  earlier  contention  that  £min  could  be  as  low  as  one  third  of  E,  due  to  the  capa¬ 
citive  structure  of  IPWs. 

We  now  consider  IPWs  launched  in  low-pressure  glow  discharges  where  the  wave  radius  is  typi¬ 
cally  restricted  to  the  glow-discharge  radius.  A  review  of  these  and  related  experiments  was  given  by 
Fowler.34  We  shall  concentrate,  however,  on  the  more  recent  experiments  of  Suzuki.19  The  latter  work 
is  selected  because  it  specifies  all  primary  parameters  of  interest.  Suzuki  additionally  presented  a 
bounded-wave  theory  similar  to  that  given  here. 

Suzuki's  results  may  be  summarized  as  follows:  wave  speed  V  scaled  as  nj/2:  V  exhibited  a 
monotonic  dependence  on  0O  which  saturated  at  high  <£n'.  v  scaled  linearly  with  where  r,  is  the 
voltage  rise-time;  and  positive  waves  typically  propagated  faster  than  negative  waves.  These  findings 
agree  well  with  the  bounded-wave  predictions.  This  agreement  was  further  confirmed  by  numerical 
simulations  which  coupled  the  chemistry  code  CHMAIR  to  Eqs.  (27)  and  (29).  using  constant  specified 
wave  speed  V and  a  constant  capacitance  of  C  -  0.2. 

9.  CONCLUSION 

In  this  paper  we  have  concentrated  on  the  propagation  of  ionizing  waves  in  preionized  gases.  We 
have  shown  that  these  waves  propagate  not  because  of  directed-particle  motion,  but  because  electron- 
impact  ionization  due  to  the  local  electric  field  raises  the  degree  of  ionization  and  hence  the  conduc¬ 
tivity  in  the  wave  head.  As  a  result  the  wave  space  charge  and  fields  are  displaced  forward,  thereby 
causing  ionization  to  commence  in  a  new  region  ahead  of  the  wave.  Processes  such  as  electron  drift, 
diffusion,  etc.  are  of  secondary  or  lesser  importance. 
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The  crux  of  the  analysis  has  been  the  use  of  similarity  and  related  arguments  to  by-pass  the 
mathematical  complexities  and  arrive  at  a  wave  velocity,  F(<£0,cr,).  The  complexities  chiefly  comprise 
two  forms:  the  two-dimensional  elliptic  character  of  Maxwell's  equations,  a  feature  which  is  responsi¬ 
ble  for  space-charge  field  enhancement  and  is  thus  essential  for  wave  propagation  in  nonuniform  fields: 
and  complications  due  to  the  functional  dependencies  of  5  and  fi  on  £,  ,V.  and  <r.  The  similarity  argu¬ 
ment  is  like  that  used  in  simple  thermodynamic  and  hydrodynamic  problems.-'5-36  It  is  based  on  the  fact 
that  the  continuity  equation  for  conductivity  o-  and  Maxwell's  equations  for  field  E  and  potential  6  are 
linear. 

In  the  simplest  nonuniform  field  problem,  we  have  found  that  wave  speed  scales  nearly  linearly 
with  the  external  control  parameters  <r,  and  4>q.  We  have  also  discussed  the  influence  of  constraints 
such  as  radial  boundaries  and  electrode  effects,  and  have  outlined  how  to  extend  the  analysis  to  propa¬ 
gation  in  un-ionized  gases.  Comparison  with  experiment  has  generally  been  favorable.  Agreement 
between  theory  and  experiment  at  a  few  operating  points  is  insufficient,  however,  to  validate  a  theory, 
as  evidenced  by  the  wealth  of  IPW  theories  and  good  agreement  claimed  for  all.  Of  greater  use  is  the 
determination  and  evaluation  of  scaling  relationships  between  measurable  variables.  Until  careful  meas¬ 
urements  such  as  those  of  Suzuki 19  have  been  performed  over  a  wide  parameter  space,  distinctions 
between  differing  theories  will  likely  remain  unresolved. 
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Appendix  A 

CHARACTERISTIC  LENGTH  SCALES 


Consider  the  following  length  scales  which  are  to  be  evaluated  at  the  peak  field.  £,: 

/,  - 

l:  - 
lx  - 

U  “ 

/<  - 

where  by  definition 

We  employ  the  reasonable  assumption  that  the  mobility  fx  and  avalanche  rate  S  also  attain  peak 
at  £,.  Hence,  at  £. 

|r  ( 5.  -  o 

and 

I  M 

In  (E/m*!  “'i- 


hi 


3  .  - 
—  !n  E 
or 

3  ,  r\ 

»7  "  C: 


~  In  <C>> 


iii  ,„r 

1  Sr  "  c= 


hi/-’ 


i_3i 

i3:: 


hi/: 


In  (/u£) 


Evaluating  Eqs.  (24)  and  (28)  at  £,  then  yields 

V  +•  /a  £,  “  S.  /j 

r  -  i,£, 

Substituting  these  results  into  Eq.  (29)  evaluated  at  £,  produces 

J_  _  J _ 1_ 

I\  h  1 1> 

where 

<b.  . 

/„-  -jr  (1  -  f°£rC.). 

Similar  manipulations  performed  on  the  first  derivatives  of  Eqs.  (24).  (28).  and  (29)  lead  to 

h  -  /. 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

( A6) 
\alues 

<A7) 

( A8) 

( A9) 
<  A 1 0 ) 

(All) 

( A 1 2 ) 

( A 1 3 ) 


19 


and 


1  J _ | _ i_ 

\+liE,/y  I ;  l i  I-' 


( A 1 4 ) 


The  dependence  of  the  length  scales  I,  on  wave  speed  Tin  Eqs.  (A11MA14)  vanishes  in  the  limit 

!/x£,l  <  <  V  <  <  ( L.C. )“,/2.  ( A 1 5 ) 

Hence,  by  imposing  this  constraint  we  can  neglect  wave  motion  and  simply  consider  the  electrostatic 
problem. 

For  a  filamentary  wave  in  which  the  radius  r„  and  body  length  /satisfy 

/•„</;<</.  ( A 1 6 ) 

as  is  implicitly  assumed  in  the  transmission-line  equations.  Gauss's  law  dictates  that  the  electrostatic 
field  near  the  tip  of  the  wave  is  determined  by  the  distribution  of  space  charge  in  the  wave  head.  Since 
the  mobility  n  is  not  a  strong  function  of  field  £.  we  thus  expect  that 

U  =  h  =  I,  <  A 1 7 » 

where 

/,  -  C.<bJE..  (A  18) 

Here  characterizes  the  space  charge  in  the  wave  head,  and  £,  is  the  peak  field  at  the  tip.  In  the 

absence  of  coaxial  ground  returns,  the  distributed  capacitance  is  given  in  Gaussian  units  by 

C  =  (2  In  (2/j/rH.)l-‘  <  A 1 9 > 

and  typically  satisfies 

C.  =  0.2.  (A  20) 


Combining  Eqs.  (A1D-1A15)  with  Eqs.  <A17).  (A18),  and  (A20)  produces  the  desired  result  that 


for  /  —  1,5. 


/,=  /,<  /  „ 


(A21) 


We  use  result  (A21)  as  follows.  Eliminating  wave  speed  I  from  Eqs.  <A9>  and  ( A 1 0 )  produces, 
to  order  /j.  EJ  K 


S  _  Ml  i  l 

'  "  C;,*,  I,  “  C,  ' 


<  A  2  2 ) 


Substituting  Eqs.  (33)  and  (34)  into  Eqs.  IA9)  and  (A22)  then  leads  directly  to  the  solution  for  wave 
speed,  T(<^,,  <t,),  as  expressed  by  Eqs.  < 35 )  and  (38). 
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Appendix  B 

ENERGY  CONSERVATION 


Previous  investigators37  -'8  have  attempted  to  derive  scaling  relationships  using  energy  conserva¬ 
tion.  We  note  here  a  potential  fallacy  in  such  an  approach,  show  how  to  avoid  it,  and  derive  an  upper 
limit  for  the  wave  velocity  that  is  more  general  and  yet  consistent  with  the  results  presented  earlier. 

The  power  supplied  by  the  electrode  goes  into  field  energy  and  gas  dynamics  according  to 

('<T  C +  T  Ml  +€.1  IBI) 

where  c,  is  the  energy  absorbed  per  unit  length  by  the  gas  and  where  the  electrode  current 

iB2) 

These  equations  reduce  to 

e,  -  1  C.6, fd  -  (B3I 

Note  that  e .  — *  0  as  l  —  <Z..C.)-1  ::  i.e..  all  the  energy  is  stored  in  the  electromagnetic  fields. 

A  certain  fraction  /  of  the  absorbed  energy  goes  into  ionization,  which  suggests  that  the  final 
wave  conductance  is  given  by 

Z,  =  /  Pfi ,  €  f/  W  (B4) 

where  IT  is  the  gas  ionization  energy.  The  wave  velocity  can  be  expressed  according  to  Eq.  (28)  by 

V  —  E,  I  ,i  <B5> 

Hence, 

r-  y  (1  -  y-L.CA.  <B6) 


The  usefulness  of  Eq.  (B6)  depends  upon  the  ability  to  estimate  /  As  we  now  show,  the  assump¬ 
tion  of  constant  /is  often  invalid. 


Integrating  the  electron  continuity  Eq.  (24)  from  the  electrode  to  the  wave  tip  yields 


£, 


*  ’  n,E, 


f  it  s  i  — — 

**  a 


(B7) 


where  we  have  ignored  the  initial  conductance.  I,,  in  comparison  with  L,.  To  relate  this  to  the  energy- 
absorbed  by  the  gas,  we  note  that  the  wave  absorbs  energy  in  the  electron  drift  frame  at  a  rate  I  E:. 
Transforming  to  the  lab  frame,  the  wave  absorbs  energy  at  a  rate 


e_  y  F,  v 

1+nE 


where  the  final  energy  absorbed  is  given  by 


I  E: 


V  +  fiE 


(B8) 


<B9> 
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We  may  thus  rewrite  Eq.  <B?>  as 


(BIO) 


I,  =  < 


_S_ 

E2 


M  / 
M 


>  e. 


where  the  average  value  of  a  variable  .v  is  defined  by 


<  x  > 


r  /J.E  bf 

) '  +  n ,  E,  X  br 


(BID 


Note  that  terms  in  parentheses  are  —1. 


Comparison  of  E qs.  (B4)  and  (BIO)  leads  to 


/  =  <^T> 

euE: 

which  re-expresses  energy  conservation.  Wave  velocity  is  thus  given  by 

E,< 2i0  <-^T>  (1  -  V2L.C.). 
nE- 

To  estimate  /  we  rewrite  Eq.  (B12)  in  terms  of  the  Townsend  coefficients  a  and  tj: 


(B12) 


(B13) 


Figure  5  shows  a  plot  of  the  fraction 


< 


a  —  t) 

eE 


lt>. 


t" 


ELZJL  ir 
eE 


(B14) 


<  B 1 5 ) 


for  weakly  ionized  air.  This  plot  assumes  an  effective  ionization  energy  of  -  15  eV,  which  equals  the 
sum  of  the  ionization  potential  plus  the  average  kinetic  energy  of  the  plasma  electrons.  The  strong 
dependence  of/'  on  field  parameter  f/.V  demonstrates  the  fallacy  of  assuming  constant  / 


Equation  ( Bf>)  can  be  used  to  prov  ide  an  upper  limit  for  wave  speed  T  by  noting  that  the  ioniza¬ 
tion  fraction 


/  <  1 . 

<Blb) 

Using  typical  values  for  the  downstream  drift  speed. 

ft,-  E,  <  10  cm/sec. 

(BIT) 

and  for  the  ionization  energy. 

W  >  10  eV. 

(BIS) 

thus  suggests  that 

l'/d>„  <  5  x  10-  cm/  Volt  -sec 

(Bid) 

irrespective  of  initial  conditions,  gas  density,  or  wave  propagation  mechanism.  This  condition  can  be 
violated  only  if  the  average  applied  field  appreciably  exceeds  the  maintenance  field  £..  or  if  the  poten¬ 
tial  wave  is  nonionizing. 


